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Abstract. Deep X-ray surveys provide the most efficient identification of 
' Active Galactic Nuclei (AGN) activity. However, recent evidence has indicated 

C^l \ that the current < 10 keV surveys do not detect the most heavily obscured 

£H . AGNs. Here we explore whether the X-ray undetected AGN population can 

be identified within the ultra-deep Spitzer survey of the GOODS-N field using 
X-ray stacking techniques. We find evidence for AGN activity in the Spitzer 

■ dataset and the strongest and hardest X-ray signal is produced by galaxies with 
£N| ' starburst-like infrared spectral slopes and median properties of z « 0.8 and 

I/ir ~ 10 11 Lq. The stacked X-ray properties suggest that obscured AGN 

■ activity is present in these sources, with a median X-ray spectral slope of T « 1 
and Lx ~ 10 42 erg s" 1 . These overall properties are consistent with the obscured 
AGN population expected to produce the unresolved X-ray background. 

. 1. Introduction 

o i 

There is a growing need for a complete census of AGN activity. The semi 
Q-f nal discovery that every massive galaxy in the local Universe harbors a super- 

q ■ massive black hole (SMBH; Mbh > 10 6 M & ) implies that all ga laxies have 

hosted AGN activity at some time during the past f» 13 Gyrs (e.g., lRees|[l984 

czj ' 
03 



Kormendv &: Richstone| [r995), To accurately trace how and when these SMBHs 



grew requires a detailed census of AGN activity that will provide, amongst other 
things, constraints on the efficiency and duty cycle of SMBH growth. The find- 
ing that the mass of the SMBH is proportional to that of the galaxy spheroid 
$_j ■ also indicates that the growth of SMBHs and their host gala xies are regulated 

in some way (e.g.. iMagorrian et al.|[l99^ ; iGebhard et al. 200dl ). Whether or not 



this growth occured concordantly is unknown but comparisons between the rate 
of star formation and SMBH growth, both in individual objects and from a 
cosmic census, will provide important insight. 



2. Selecting and Identifying AGN Activity 

Since obscured AGNs outnumber unobscured AGNs by a factor of ~ 3-10, any 
"complete" census of AGN activity must be sensitive to the detection and iden- 
tification of obscured sources. X-ray emission is relatively impervious to obscu- 
ration, at least at high energies, and the detection of luminous or hard X-ray 
emission provides an almost unambiguous identification of AGN activity even 
in the absence of source redshifts. For this reason, the most complete census of 
AGN activity to date is provided by deep X-ray surveys. Fo r example, the 2 Ms 
Chandra Deep Field North (CDF-N; lAlexander et al.il2003h survey identifies a 
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Figure 1. (a) observed X-ray luminosity vs redshift for the spectroscopi- 
cally identified X-ray sources in the CDF-N/GOODS-N field. The tracks of 
the Compton-thick but intrinsically luminous (Lx ~ 10 44 erg s _1 ) AGNs in 
NGC 6240 (long dash) and NGC 1068 (short dash) are shown to indicate their 
detectability at different redshifts. The shaded regions show where the spec- 
troscopically unidentified X-ray sources are likely to lie ; ss 50% of the X-ray 
AGNs (almost all at z > 1; e.g., lAlexander et all [20011 : iMainieri et al.ll2005f ) 
still lack spectroscopic redshifts. (b) rest-frame far-IR luminosity vs redshift 
for spectroscopically identified 24 /im Spitzer sources in the GOODS-N field; 
the far-IR luminosities are calculated from the 24 /im data using the average 
of four SED templates (M 51, M 82, Arp 220, Mkn 231). Note that while 
NGC 1068 will only be detected out to z « 0.4 in the X-ray observations, it 
should be detectable out toz« 1-2 in the Spitzer observations. 



> 10 times larger AGN source density (« 7200 dee" 2 . iBauer et all 12001 ) than 
found in optical surveys; see Fig. la for the sensitivity of the CDF-N to AGN 
objects. 

However, although deep X-ray surveys have proved st unningly effective at 
ident ifying both obscured and unobscured AGN activity (see lBrandt fc Hasinger 
120051 for an overall review), there is clear evidence that a large fraction of the 
AGN population remains undetected in the < 10 keV (observed frame) band 
probed b y these surveys : (1) a bout half of the X-ray background is unresolved at 

> 6 keV ( Worslev et al.l 12005). (2) the obscured: unobscured AGN ratio is much 
lower than that found for comparably luminous AGNs in the local Universe (e.g., 
Treister fc Urrvl 120051 ) . and (3) few Compton-thick AGNs (N u > 10 24 cm" 2 ) 
have been identified, even though they comprise ~ 50% of the AGN population 
in the local Universe (e.g., iTozzi et al J 12009 : iPaee et alJl2003 ). 

Many of these X-ray undetected AGNs are likely to be intrinsically lumi- 
nous sources that are heavily obscured (iVn > 3-10 x 1 23 cm" 2 ) by an opt ically 
and geometrically thick torus of gas and dust (e.g., I Worslev et al.l 120051 : see 
Fig. la for the detectability of the intrinsically luminous AGN in NGC 1068 
in the CDF-N observations). Although the primary emission from these ob- 
scured AGNs will be absorbed by the torus and rendered weak or invisible at 
X-ray-optical wavelengths, it will heat the dust grains and re-emit this emis- 
sion relatively isotropically in the infrared (IR) waveband (rest-frame > 2 ^m; 



e.g.. iGranato Danesd[T994l ; lEfstathiou <fe Rowan-Robinsonlll995l ). The intrin- 



sically luminous but X-ray undetected AGNs should therefore be detected in 
deep IR surveys; see Fig. lbQ Since star formation can also produce luminous 
IR emission it is necessary to distinguish between AGNs and starburst galaxies. 
A number of recent studies have demonstrated that IR colours can often distin- 
guish between AGN and star-formation activity using crude IR dat a, with the 



form e r having "hotter" IR spectral energy distributions (SEDs; e.g.. lLacv et al.l 



20041 : IStern et alJl2005l ). However, since a large fraction of the AGN popula- 



tion in the local Universe has IR colors similar to those of starburst galaxies, 
other indicators (e.g., the detection of high-excitation emission lines, hard X-ray 
emission from stacking analyses of samples of objects, and luminous unresolved 
radio emission) are also required to identify AGN activity from star formation. 

The primary aim of this article is to investigate the prospects for select- 
ing X-ray undetected AGN with Spitzer observations in the ultra-deep CDF- 
N/GOODS-N field to increase the AGN census beyond that identified from the 
X-ray data alone. Here we select X-ray undetected AGN candidates using the 
Spitzer GOODS-N observations and then search for the presence of hard X-ray 
emission using X-ray stacking techniques. From these results we place con- 
straints on the demographics and properties of AGN activity in the IR galaxy 
population. 



3. Searching for a large X-ray-undetected AGN population in the 
CDF-N/GOODS-N field 

The Spitzer GOODS observations (P.I. M. Dickinson) were performed as part 
of the Spitzer Legacy program and cover ~ 160 arcmin 2 of the most sensitive 
regions in each of the Chandra deep fields (the CDF-N and the CDF-S). The 
IRAC 3.6-8 /im and MIPS 24 //m Spitzer observations are the deepest across 
the sky, achieving typical 1 a sensitivities of w 0.03-0.3 //Jy (3.6-8 /im) and 
« 5 fiJy (24 /mi). Here we focus on the CDF-N/GOODS-N field, which has 
the deepest X-ray coverage to date. To provide maximum sensitivity to X-ray 
faint AGN activity we restrict our analyses to the region within 5.5 arcmin of 
the Chandra aimpoint (95 arcmin 2 ), wher e the X-ray sens it ivity and angular 
resolution is optimal; see Figs. 18 & 19 of lAlexander et ah ( 20031 ) . There are 



223 X-ray sources, ~ 1800 MIPS 24 fxm sources, and ~ 3400 IRAC 8 /im sources 
within this region. 

3.1. X-ray Detected Infrared Galaxies 

The X-ray sources were matched to t he IR sources on the b asis of the X-ray 
source positional uncertainties given in Alexander et al. ( 20031 ). The positional 



1 High-resolution radio observations also provide effective obscuration-independent selection of 
AGNs. However, most radio surveys are either confined to relatively rare AGN source pop- 
ulations, such as high-z radio-loud AGNs or to AGN activity in nearby galaxies (but see 

iDonlev et alJl2005l ). Future radio observatories such as the Square Kilometer Array will pro- 
vide an extremely sensit ive AGN census but will not be fully operational for > 10 years (e.g., 

Uarvis fc Rawling3l2004f ); see |http://w ww.skatelescope.org"7| 
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Figure 2. 24 /im flux density vs 0.5-8.0 keV flux (top), X-ray detected 
AGN fraction (middle), and X-ray detected starburst fractio n (bottom). The 
shade d region shows the flux ratios of X-ray bright AGN (|Piccinotti et al.l 
119821 ) extrapolated to fainter fluxes. The decrease in the X-ray detected AGN 
fraction towards fainter 24 /im flux densities is at least partially due to a 
relative decrease in the X-ray-IR flux ratio (i.e., there is less area below the 
shaded region at faint 24 /im flux densities) and suggests that the majority 
of the X-ray undetected AGN population will have faint IR fluxes. 



uncertainites were linearly increased to give a 3% spurious source match frac- 
tion; extensive Monte-Carlo tests were used to determine the spurious match 
fraction. Overall, 78% and 87% of the X-ray sources were matched to 24 fj,m 
an d 8 /im sources, res pectively. Using the multi-wavelength source classification 
of lBauer et all (|2004h we distinguished between X-ray detected AGNs and X- 
ray detected starbursts. The X-ray sources comprise 9.4% of the 24 /im source 
population (7.0% X-ray AGNs and 2.4% X-ray starbursts) and 5.6% of the 8 /im 
source population (4.3% X-ray AGNs and 1.3% X-ray starbursts). 

Fig. 2 shows the X-ray vs 24 /im properties of the X-ray sources, including 
the X-ray AGNs and X-ray starburst fraction for different 24 /im flux density 
bins. As expected, the majority of the X-ray AGNs have X-ray-IR flux ratios 
higher than those of the X-ray starbursts. Fig. 2 also shows a clear drop in the 
X-ray detected fraction of IR sources towards the faintest IR fluxes. This is at 
least partially due to the linear decrease in the observed X-ray-IR flux ratio 
with 24 /im flux density, suggesting that the majority of the X-ray undetected 
AGN population will have faint IR fluxes0 



2 This is clearly a simplistic prediction since Fig. 2 doesn't take into account possible changes 
in redshift and SED for sources in the different 24 (im flux-density bins but it is at least 
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3.2. Searching for X-ray Undetected AGNs 

Evidence for an X-ray undetected AGN population can be explored using X-ray 
stacking techniques. The principle behind this approach is to stack the X- 
ray data of individually undetected Spitzer sources and determine the probabil- 
ity that the stacke d X-ray signal could be produced by the background (e.g., 
Brandt et al.|[200ll ; the background is determined from a Monte-Carlo analysis 



of source-excluded regions around the stacked sources. The expected signature 
of a significant obscured AGN population is a detection in the hard (2-8 keV) 
band and a comparatively flat X-ray spectral slope (r < 1.4). By comparison, 
unobscured/low-luminosity AG N and typical starburst galaxies hav e softer X- 
ray spectral slopes (T « 2; e.g.. iPtak etaDll999l : [George et alJ^OOOT ). We note 
that neutro n-star X-ray binary populations can also have flat X-ray spectral 
slopes (e.g., IColbert et al.l [2004 ) . however, it seems unlikely that the integrated 



X-ray emission from a starburst galaxy will be dominated by this comparatively 
rare galactic population. 

For the stacking analyses we have used the procedure outlined in lWorslev et al 



(2006) and stacked Spitzer galaxies in the soft (0.5-2 keV), hard (2-8 keV), and 



full (0.5-8 keV) bands. The background was determined using 1,000 Monte- 
Carlo trials for each stacked source and the X-ray counts were measured using 
apertures based on the 90% encircled-energy radius. Stacking all of the X-ray 
undetected IR sources we obtain highly significant detections (> 10 er) in all 
three X-ray bands, corresponding to an overall X-ray spectral slope of T = 1.4 
and suggesting that some of the stacked sourc es are likely to be obsc ured AGNs 
(contrast with the T = 2.0 result obtained bv lAlexander et al.ll2002l using shal- 
lower 15 /im data over a smaller area); this is further confirmed by a 3.5 a 
detection in the ultra-hard (4-8 keV) band. We also stacked the data in the 
6-8 keV band and obtained a non detection in this band; the 3 a upper limit 
corresponds to < 20% of the 6-8 keV background. If a substantial fraction of the 
X-ray background is unresolve d in the 6-8 keV band by current X-ray surveys, 
as suggested by the analyses of lWorslev et al.l (|2005l ). then the other contribut- 



ing sources must be either IR-faint (/24^m < 25 /^Jy) or comparatively rare in 
~ 100 arcmin 2 regions. 

A basic interpretation of the data in Fig. 2 is that X-ray undetected AGNs 
will be predominantly detected at the faintest IR fluxes. To test this simple pre- 
diction we stacked the X-ray data for individually undetected sources in different 
24 fim flux-density bins; see Fig. 3a. The 2-8 keV stacked flux generally increases 
towards fainter 24 /iin flux densities, suggesting an increase in the contribution 
from obscured AGN activity, with the /24^ m = 75-225 [iJy bin producing the 
strongest and hardest X-ray signal. However, the softer X-ray spectral slope for 
sources with f^Aum = 25-75 //Jy implies that starburst galaxies probably domi- 
nate the stacked X-ray emission at the faintest IR fluxes. Furthermore, the com- 
paratively weak 2-8 keV flux suggests that there may be fewer obscured AGNs 
than at fiA^m = 75-225 /iJy, possibly indicating a genuine drop in the number 
of obscured AGNs at low IR luminosities or at high redshifts {f2A^m > 75 fiJy 



qualitatively s i milar to the conclusions drawn from the AGN population synthesis model of 
iTreister et all (|2006D . 
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Figure 3. (a) stacked 2-8 keV flux density vs 24 /im flux density for indi- 
vidually undetected Spitzer galaxies; the X-ray spectral slope for each 24 /im 
flux-density bin is also shown. The strongest and hardest X-ray signal is 
produced by /24/xm = 75-225 /iJy sources. The dashed line indicates the over- 
all stacked X-ray signal obtained by stacking all of the Spitzer galaxies, (b) 
stacked 2-8 keV flux density vs X-ray spectral slope for individually unde- 
tected Spitzer galaxies with different IR spectral slopes (a) . The results for 
the a = 0.5-1.0 sources are also given for different flux ranges (/24pm = 25- 
75 /iJy and f^ium > 75 /iJy); the strongest and hardest X-ray signal comes 
from the f^Ajim > 75 /LtJy sources. 



corresponds to Lfir > 10 11 L at z = 1). These results are qualit atively similar 
to tho se expected from the AGN population synthesis model of lTreister et al ' 



2006) and the interpretation of shallower Spitzer and X-ray data bv lBrand et al 



20061 ). More detailed analyses are required to provide definitive conclusions. 



We can also explore the average X-ray properties of the Spitzer sources as 
a function of IR spectral slope (a, where f u a v a ). In the local Universe many 
AGNs have flat IR spectral slo pes while starburst galaxies a re often distinguished 
by steeper IR spectral slopes. lAlonso-Herrero et al.l (|2006l ) investigated the de- 
tailed properties of Spitzer galaxies with flat IR spectral slopes (a < —0.5) in 
the 1 Ms CDF-S field and deduced that only « 50% of these IR-selected AGNs 
are detected in the X-ray band. Although the flat IR spectral slopes of these 
sources clearly indicate that they host AGN activity, they are comparatively 
rare (« 800 sources deg -2 , including both X-ray detected and X-ray detected 
AGNs) and unlikely to dominate the unresolved > 8 keV background. As a 
comparison, the range of IR spectral slopes for the X-ray detected AGN in the 
CDF-N/GOODS-N dataset is a = -2.0-1.0, and only « 30% have a < -0.5, 
indicating that many have IR spectral properties more consistent with those 
of starburst galaxies. Since the aim of the study here is to provide the most 
complete census of AGN activity, we have stacked the X-ray data of Spitzer 
galaxies over a wide range of IR spectral slopes (a < —0.5, a = —0.5-0.0, 
a = 0.0-0.5, a = 0.5-1.0). The X-ray spectral slope and 2-8 keV flux density 
of the stacked sources are shown in Fig. 3b. All are significantly detected in 
the 2-8 keV band but the strongest and hardest X-ray signal comes from the 
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sources with the steepest IR spectral slope (a = 0.5-1.0). The steep IR spectral 
slope of these sources does not obviously imply the presence of AGN activity 
but the hard X-ray signal (T = 1.2) indicates that a considerable fraction of the 
sources host obscured AGN activity. Motivated by the discovery that sources 
with f2i^m > 75 /zJy show greater evidence for X-ray undetected AGN activity 
(see Fig. 3a), we further stacked the sources in two different 24 [im flux bins 
(/24/xm = 25-75 /iJy and /24/xm > 75 /uJy). The flatter X-ray spectral slope 
(r = 1) and larger 2-8 keV flux density of the fiA^m > 75 /iJy sources with 
a = 0.5-1.0 clearly indicates that the majority of the obscured AGNs reside in 
this population. 



4. Discussion 



The X-ray stacking analyses provide evidence that a large X-ray undetected 
AGN population does indeed reside in the Spitzer galaxy population. The domi- 
nant AGN population has an IR spectral slope of a = 0.5-1.0 and fv&um > 75 /iJy. 
Their IR colors are more similar to starburst galaxies than AGNs but X-ray 
stacking analyses indicate that they host obscured AGN activity. We can gain 
further insight into the properties of these sources by investigating their median 
redshifts, luminosities, and source densities. 

The median spectroscopic redshift of the a = 0.5-1.0 Spitzer galaxies is 
z ~ 0.8, which corresponds to an median luminosity of Lftr ~ 10 11 L & ; the 
majority of the redshifts were taken from Wirth et al. ( 2004! ) . Assuming this 



median redshift, the 2-8 keV luminosity of the stacked sources with a = 0.5-1.0 
and /24^tm > 75 fiJy is ~ 10 41 erg s _1 (rest-frame 3.6-14.4 keV). This X-ray 
luminosity should be considered a lower limit since (1) many of the stacked 
galaxies will not be AGNs and their inclusion in the stacking analysis will de- 
crease the average 2-8 keV flux, and (2) the luminosity has not been corrected 
for the presence of instrinsic absorption. These two issues are briefly considered 
below. 

We cannot directly determine the "dilution" to the hard-band X-ray sig- 
nal from non-AGN activity, however, we can place constraints from the non- 
detection of individual sources in the X-ray data. For example, if the X-ray 
undetected AGN represented just 1% of the stacked IR galaxies but contributed 
all of the 2-8 keV flux then they would have individual 2-8 keV luminosities 
~ 2 orders of magnitude higher and would be individually detected in the X-ray 
data. Therefore, since the stacked 2-8 keV flux is an order of magnitude below 
the individual X-ray sensitivity limit, these X-ray undetected AGNs are likely to 
represent > 10% of the stacked IR galaxies. This indicates an average observed 
(i.e., uncorrected for absorption) 2-8 keV luminosity of ~ 10 41 -10 42 erg s _1 and 
places lower limits on the X-ray undetected obscured AGN source density of 
» 500-5,000 deg~ 2 (for AGNs with a = 0.5-1.0 and / 2 4 M ™ > 75 /iJy). 

A crude constraint on the amount of absorption towards the AGNs can 
be determined from the average X-ray spectral slope for a given AGN model 
(e.g., an observed X-ray spectral slope of T = 1.0 from a z ~ 0.8 AGN would 
imply Ah ~ 3 x 10 22 cm s~ 2 for a simple intrinsic AGN power-law spectrum 
of r = 2.0). Since the soft-band flux almost certainly has a large contribution 
from star-formation activity, this approach effectively provides a lower limit on 
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the true absorbing column density. Another approach, albeit very uncertain, is 
to estimate the intrinsic (i.e., corrected for absorption) X-ray luminosity of the 
AGNs from the IR luminosities and then compare this to the observed X-ray 
luminosity Under the assumption that 10-100% of the IR luminosity is powered 
by an AGN , the median X-ray luminosity of these sources estimated from the 
lElvis et ail ( 19941 ) SED is « 10 42 -10 43 erg s _1 , suggesting large amounts of 
absorption (JV H » 10 23 cm -2 ). 

Although these analyses on the properties of the X-ray undetected AGNs 
are crude, their estimate d prop erties are qualitatively consistent with those re- 
quired by IWorslev et al.1 ( 20051 ) to produce the unresolved X-ray background 
(i.e., z w 0.5-1.5, N H « 10 23 -10 24 cm" 2 , L x « 10 42 -10 43 erg s" 1 ). Interestingly, 
the median redshift, IR luminosity, and X-ray luminosity range of the X-ray de- 
tected AGNs with a = 0.5-1.0 are also similar to those estimated for the X-ray 
undetected AGNs with a = 0.5-1.0. This is entirely consistent with that ex- 
pected for the unified AGN model, where the X-ray undetected AGNs represent 
the more heavily obscured counterparts of the X-ray detected AGNs. Detailed 
observations of individual objects are now required to elucidate the properties 
of this X-ray undetected AGN population. 
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